Tunable Fermi acceleration in the driven elliptical billiard 
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Abstract 

We explore the dynamical evolution of an ensemble of non-interacting particles propagating freely 
in an elliptical billiard with harmonically driven boundaries. The existence of Fermi acceleration 
is shown thereby refuting the established assumption that smoothly driven billiards whose static 
counterparts are integrable do not exhibit acceleration dynamics. The underlying mechanism 
based on intermittent phases of laminar and stochastic behavior of the strongly correlated angular 
momentum and velocity motion is identified and studied with varying parameters. The diffusion 
process in velocity space is shown to be anomalous and we find that the corresponding characteristic 
exponent depends monotonically on the breathing amplitude of the billiard boundaries. Thus it is 
possible to tune the acceleration law in a straightforwardly controllable manner. 

PACS numbers: 05.45.-a,05.45.Ac,05.45.Pq 
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Fermi acceleration was first proposed in 1949 |l| to explain the high energies of cosmic ray 
particles interacting with a time- dependent magnetic field (for a review see Nowadays, 
Fermi acceleration is investigated in a variety of systems belonging to different areas of 



physics, such as astrophysics 



|9| and has even 



m 

, plasma physics , atom optics 
been used for the interpretation of experimental results in atomic physics 10|. In the context 
of dynamical system theory, Fermi acceleration is defined as the unlimited growth of the 
energy of a particle moving in a time- dependent potential. For the experimental realization 
of Fermi acceleration in the laboratory, driven, i.e. time-dependent billiards are the most 
promising devices. A particle propagating in these systems gains or looses energy whenever 
it hits the billiards' moving boundary. An infinite sequence of such collisional events can 
in principle lead to Fermi acceleration of the scattered particle. Driven billiards can be 
realized experimentally in various ways, such as (superconducting) mesoscopic cavities (see 
Il2| and refs. therein) or atom optical setups employing acousto-optical scanners 



ref. 

for the (time-dependent) deflection of correspondingly detuned laser beams |13[ |. 



The simplest realization of a driven billiard is the well-known Fermi-Ulam model (FUM) 
[isl . 16 1 which consists of an ensemble of non-interacting freely moving particles in one 



dimension bouncing between an oscillating and a fixed wal 
have been the subject of extensive theoretical (see Re: 



experimental 



13, 



18 



19| studies. It has been proven 



The FUM and its variants 



15l | and references therein) and 



16l | that a necessary condition for 



the development of Fermi acceleration in the FUM is the non-smooth (non-differentiable) 
dependence of the velocity of the moving wall on time (this includes especially randomized 
dynamical systems). An alternative scenario for the occurrence of Fermi acceleration in ID 
systems is a bouncer model 20| with an additional nonlinear time-independent term in the 
potential. 

In contrast to the ID case, driven billiards with a higher dimensional configuration space 
pose many open questions. The lack of the possibility to visualize phase space in terms of 
^oincare surfaces of section complicates the analysis of e.g. 2D time-dependent billiards 



2l| . Despite this fact, several studies of Fermi acceleration in these systems 22, 



3, 



24 



25| 



have been performed. These investigations provide us with the conclusion that a sufficient 
condition for the occurrence of Fermi acceleration in a 2D smoothly driven billiard is the 
existence of a chaotic part in the phase space of the corresponding time-independent system, 
obtained by assuming static boundaries (see the LRA conjecture in ref. 23|). This is 
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supported by the absence of Fermi acceleration in the smoothly oscillating circular billiard 
where the corresponding static system is integrable |26| . On the other hand the static Lorentz 
gas possesses a predominantly chaotic phase space and its harmonically driven counterpart 



exhibits Fermi acceleration 



24 



271] . Similar results were obtained for the oscillating stadium- 



like oval 291] as well as annular billiard 25|], each of them possessing chaotic portions 
with respect to the phase space of the corresponding time-independent system. 

To elucidate the above-described problem of the existence of Fermi acceleration for driven 
billiards whose time-independent counterparts are integrable, we investigate here the ellip- 
tical harmonically driven billiard. According to the existing literature the ellipse should not 
exhibit Fermi acceleration (see e.g. ref. 2^). We show that, opposite to what is expected in 
the literature, Fermi acceleration occurs for ensembles of particles propagating in a breath- 
ing elliptical billiard even for very small oscillation amplitudes. The driving causes a layer of 
unstable motion around the separatrix of the elliptical billiard leading to large fluctuations 
of the velocities of the particles as they cross it. These fluctuations, which increase with 
time, lead to an anomalous diffusion in velocity space and an unlimited growth of the mean 
kinetic energy of the multiple scattered particles. The anomalous behavior of the velocity 
diffusion process is due to the intermittent character of the dynamics attributed to the in- 
stability of the fixed points in the regime of librators. We determine the dependence of the 
exponent of the acceleration law on the amplitude of the oscillation of the boundary leading 
to the conclusion that the smoothly driven elliptical billiard represents a tunable source of 
Fermi acceleration. 

Let us specify our setup (for a more detailed description see ref. |3l|)- The static ellipse 
is integrable: besides the energy E the product of the angular momenta F about the two 



foci of the ellipse 



30| is conserved. 



p2(l + (l-£2)cotV)-£' 



(1) 



1 + (1 - £2) cot2 - £2 ' 

where e is the numerical eccentricity. The variables p = cos a and are defined at the 
elliptical boundary of the billiard: a is the angle between the tangent and the trajectory of 
the colliding particle while Lp is the azimuthal angle, see eq. ([2]) and Fig. [H The phase space 
is globally divided by the separatrix {F = 0), associated with two hyperbolic fixed points, 
into rotators {F > 0) and librators F < 0. Two elliptic fixed points are located at the 
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global minimum Fmin of F. The driving law of the elliptic boundary is: 

x{t) I I {Aq + C sm{ujt)) cos (f I 
y{t) I I (5o + Csin(a;t)) sinv9 I 

where t is time, {x{t),y(t)) is a point on the boundary, ip is a 27r-periodic parameter, C > 
is the driving amplitude, u is the frequency of oscillation and Aq, Bq are the equilibrium 
values of the long and the short half-diameter, respectively. The phase of oscillation at t = 
is set to zero in the following. In order to reduce the number of parameters to be varied 
in our investigations, we fix a; = 1, = 2 and Bq = 1 (arbitrary units). The scattering 
dynamics in the billiard is described by an implicit AD map specifying the sequence of the 
coUisional events on the boundary of the ellipse. A convenient choice of the variables of this 
map is {Lpn,pn, \vn\,^n), wherc \vn\ is the magnitude of particles' velocity after the collision 
and C,n is the phase of the boundary oscillation at the n-th collision. The implicit equations 
defining this map are: 

<(t„+i - t„) + / vy{tn+i - t„) + y,^'^ 



_^ -nV-^+i -^nj I yn _ 1 = Q, (3a) 

Aq + C sin(cjt„+i) J \Bo + C sin(a;t„,+i) / 

Xn+l = Xn + Vn{tn+1 ~ tn) (3b) 
Vn+1 =Vn-2 [hn+1 " (v„, - Un+l)] ' n„+i (3c) 

where the smallest tn+i > tn that solves (ISal) has to be taken and ^„,+i = tn+i mod 27r. 
Xn = {xn,yn) IS the nth colhsion point, (pn+i can be obtained by inverting ([2]). In eq. 
(!3cl) . Un+i is the boundary velocity and n„+i the normal vector of the collisional event 
occurring at time tn+i and position Xn+i. When iterating the mapping ([3]) numerically, 
solving (|3al) requires the major computational effort in terms of CPU-time, even when 
applying indispensable advanced bracketing techniques. When the driving is applied, both 
E and F are no longer conserved quantities. However, they are still very valuable for the 
description of the dynamics in the presence of the driving ([2]), since at every time instant 
the elliptical form of the boundary is preserved. 

Before discussing the ensemble averaged properties of the time-dependent system, it is 
illuminating to analyze a typical trajectory being initially on a librator orbit (C = 0.2). 
Fig. [2] shows the evolution of F (upper curve, light gray) and \v\ (lower curve, light gray) 
respectively as a function of the number of collisions n. According to Fig. [2^, F{n) alternates 
between periods of regular (laminar phases) oscillations (intervals [0, 2.8 ■ 10^] and [1.4 • 
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FIG. 1: Schematic picture of the setup of the ehiptical bihiard. 

10^, 1.75-10^]) and periods of irregular fluctuations (turbulent phases). During these laminar 
phases, F{n) never crosses the F = line, whereas during the turbulent phases F{n) remains 
essentially within the zone [Fmin,Fmax] = [—1,0.4], repeatedly crossing the F = line 
associated with the separatrix (the separatrix is defined by F = 0, independent of t). From 
Fig. [2]d and [2t we see that this structure of laminar and turbulent phases exists on different 
scales of n. Whenever F{n) is in a laminar phase, I^TI oscillates around a fixed central value, 
whereas during the turbulent F{n) periods, \v\{n) starts to develop intervals with strong, 
irregular fluctuations leading to a sudden increase or decrease of its central value. In fact, 
the increasing parts of the trajectory |'y|(?T.) prevail such that a net increase of the velocity 
for longer times can be observed. 

The above-observed behavior indicates two important dynamical properties of the driven 
system: (i) the region around the separatrix is characterized by stochastic dynamics and 
(ii) the trajectories F{n) and \v\{n) are strongly correlated (see Fig. [2]). The enhanced 
stochasticity of the dynamics around F is also verified by the appearance of the corre- 
sponding power-spectrum S{k) = ^ J2n=i F{n)e '^~n^ ^ presented in Fig. O Clearly, 
the continuous part of the spectrum is dominating in this case (Fig. [8)3). On the contrary, 
in regions of large negative F values the power spectrum, shown in Fig. [3^, is dominated 
by isolated peaks while the continuous background is much less pronounced. Furthermore, 
the correlated behavior of |i7| and F suggests that an acceleration mechanism may be as- 
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sociated with the crossing of the F — line. To elucidate this, it is necessary to integrate 




FIG. 2: The functions F{n) (upper curve, light gray) and \v\{n) (lower curve, light gray) for a 
single trajectory for C = 0.2 in arbitrary units. The solid black lines indicate the cumulative mean 
in each case. The two successive magnifications demonstrate that the corresponding fluctuations 
occur on many scales. 

out regular as well as stochastic fluctuations of the laminar phase occurring on small time 
scales. For such an averaging procedure to be successful one has to detect first efficiently the 
intervals of laminar evolution. A rough estimation of these intervals is obtained by analyzing 
the timeseries F{n). As 'working intervals' (to calculate the cumulative mean) resembling 
the laminar phases, we identify the dynamics between two successive zero crossings of F. 
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Then, in each such interval j, we calculate the cumulative mean F{n)j = ^ ^2^=1 

— n^^=i ^ function of n. Here n takes values n = 1, ..,Lj with Lj being 

the number of collisions between the j — 1-th and the j-th zero crossing of F. The results 
for Fj and \v\j of the above-considered trajectory are shown as solid black lines in Fig. [21 
The similarity of the dynamics of the cumulative mean of F with a typical trajectory show- 



ing intermittency 



32| is remarkable. Additionally, the one-to-one correspondence between 



periods of stochastic behavior (bursts) of Fj and periods of fluctuations of \v\j is obvious. 
The intermittent character of the F-dynamics goes along with a specific distribution of the 
laminar lengths Lj, i.e. the intervals between successive zeros of F, which according to 



the theory of intermittent maps 



331] should obey a power-law. Indeed, we find asymptoti- 



cally (for large values of Lj) a power-law behavior of the distribution with a characteristic 
exponent of = —2.00 ± 0.05. 



(7) 



10^ 



10^ 



10^ 



10" 




CO 



10^ 



10^ 



10^ 



10" 



0.6 0.8 




FIG. 3: Power spectrum S{k) of a typical laminar (a) and turbulent phase (b) in arbitrary units. 

Let us now investigate whether Fermi acceleration occurs in the driven ellipse by exam- 
ining {\v\){n), where ( . ) denotes ensemble averaging. In Fig. H] we show the results of 
our simulations obtained using 1000 trajectories with initial values of a and Lp uniformly 
distributed in [0,7r] and [0,27r], respectively, for different driving amplitudes C (see also Ta- 
ble I). Each trajectory is propagated 10^ collisions and possesses the initial velocity \vq\ = 1. 
Obviously, we encounter Fermi acceleration for all ensembles, i.e. driving amplitudes consid- 
ered here. After an initial transient behavior of typically O(IO^) collisions characterized by 
relatively small variations of the mean velocity, the function {\v\){n) becomes an increasing 
power-law {\v\)(n) ~ n^^*^) {P{C) > 0). The associated exponent (3, as displayed in Table 
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FIG. 4: The function (|'C|)(n) (arbitrary units) for an ensemble of 10'^ trajectories for different 
values of the driving amplitude C, together with the best linear fit (between the two crosses). 



TABLE I: Exponent /3(C) for (|'w|)(n) ~ n^(<^) for different values of the amplitude C. 
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I, depends monotonically on C. The diffusion in velocity space is anomalous in agreement 
with the fact that the dynamics in F-space is intermittent. This statement holds, of course, 
not only for (|'u|)(n) but also for the dependence on the real time (|'u|)(t). The above analysis 
allows us to identify the basic ingredients of the mechanism responsible for Fermi accelera- 
tion in the driven elliptical billiard. For C = the voluminous librator region in phase space 
contains a continuous set of invariant and impenetrable manifolds. In the presence of the 
driving (C ^ 0) it becomes penetrable and individual trajectories starting in the librator 
part of phase space are 'pushed' towards the separatrix (ejection phase). Due to the driving 
the neighborhood of the separatrix is replaced by a chaotic layer whose properties depend 
on the amplitude C of the boundary oscillation. While the dynamics is laminar in the ejec- 
tion phase, the motion in the chaotic layer is to a large extent stochastic, introducing large 
fluctuations of F and consequently of the velocity which in turn leads to Fermi acceleration. 
In conclusion we have investigated the development of Fermi acceleration (FA) in the 
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harmonically driven elliptical billiard. Opposite to the expectations according to the state- 
of-the-art of the field, we discover FA of the form {\v\){n) ~ with the exponent (] being 
tunable via the amplitude of the oscillation of the elliptic boundary. We remark that no FA 

□ n 

has been observed in Ref. [22| in the driven ellipse. However, Ref. [22] employs a slightly 
different driving mode. Additionally, the transient in which there is no FA is rather long and 
could be easily misinterpreted as the absence of FA. Here we demonstrate for the first time 
that FA is observed in a higher dimensional driven system for which the corresponding static 
counterpart is integrable. To our knowledge, the considered system is also the first billiard 
system showing an acceleration law that is tunable in a straightforward and controllable 
manner despite the complex structure of the underlying higher dimensional phase space. 
We show that the observed tunable anomalous diffusion in velocity space is associated with 
an intermittent dynamics of the second integral (F) of the static integrable system, all by 
all triggered by the driving. 
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